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1 Introduction 



The study of the fundamental dynamical processes associated with the propagation of 
high-power ultrashort pulses in optical fibres is of paramount importance: the non-linear 
soliton(ic), or near-soliton(ic), operating mode(s) in such systems are very promising; in 
particular, very high data transmission rates, high noise immunity, and the accessibility of 
new frequency bands Q. The classical, mathematical model for non-linear pulse propaga- 
tion in the picosecond time scale in the anomalous dispersion regime in an isotropic, ho- 
mogeneous, lossless, non- amplifying, polarisation-preserving single-mode optical fibre is the 
non-linear Schrodinger equation, NLSE |fl, 0; however, in the subpicosecond-femtosecond 
time scale, experiments and theories on the propagation of high-power ultrashort pulses 
in long monomode optical fibres have shown that the NLSE is no longer valid and that 
additional non-linear terms (dispersive and dissipative) and higher-order linear dispersion 
should be taken into account [||. In this case, subpicosecond-femtosecond pulse propagation 
is described (in dimensionless and normalised form) by the following non-linear evolution 
equation (NLEE) [§, 

1 ~ ~ T 

id^u + -d 2 u + \u\ 2 u + isd T (\u\ 2 u) = —iTu + i5d^u + — ud T (\u\ 2 ), (1) 

2 To 

where u is the slowly varying amplitude of the complex field envelope, £ is the propagation 
distance along the fibre length, r is the time measured in a frame of reference moving with 
the pulse at the group velocity (the retarded frame), s (> 0) governs the effects due to the 
intensity dependence of the group velocity (self-steepening), T is the intrinsic fibre loss, 5 
governs the effects of the third-order linear dispersion, and t„/to, where tq is the normalised 
input pulsewidth and r n is related to the slope of the Raman gain curve (assumed to vary 
linearly in the vicinity of the mean carrier frequency, ujq), governs the soliton self- frequency 
shift (SSFS) effect g. 

Our strategy to the study of (1) will be the following: (i) by setting the right-hand 
side of (1) equal to zero, we obtain the following equation (integrable in the sense of the 
inverse scattering method (ISM) 

id t u + ^d 2 u + \u\ 2 u + isd x (\u\ 2 u) = 0, (2) 

which we hereafter call the modified non-linear Schrodinger equation, MNLSE (note that, 
the physical variables, £ and r, have been mapped isomorphically onto the mathematical t 
and x variables, which are standard in the ISM context); and (ii) since T, 5, and t u /tq are 
small parameters Q, we treat (1) as a non-integrable perturbation of the MNLSE. 

From the above discussion, it is clear that (multi-) soliton solutions of the MNLSE 
play a pivotal role in the physical context related to optical fibres. However, since prac- 
tical lasers cannot be designed to excite only the soliton(ic) mode(s), but also excite an 
entire continuum of linear-like dispersive (radiative) waves, to have physically meaningful 
and practically representative results, it is necessary to investigate solutions of the MNLSE 
under general initial (launching, in the optical fibre literature [|J) conditions, without any 
artificial restrictions and/or constraints, which have both soliton(ic) and non-soliton(ic) 
(continuum) components; physically, it is towards such a solution that the initial pulse 
launched into an optical fibre is evolving asymptotically ||. Also, one may think of general 
initial conditions (Cauchy data) as a perturbation of the so-called reflectionless potential 
(soliton(ic)) initial conditions, and, since the action of any perturbation is always more or 
less non-adiabatic, the soliton(ic) component of the solution will be accompanied by the 
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appearance of a non-soliton(ic) component [fj]. From the physical point of view, therefore, 
it is seminal to understand how the continuum and the soliton(s) interact. For several 
soliton-bearing equations, e.g., KdV, Landau-Lifshitz, and NLS, and the RMB (reduced 
Maxwell-Bloch) system, it is well known that the dominant (0(1)) asymptotic (i^oo) ef- 
fect of the continuous spectra on the multi-soliton solutions is a shift in phase and position 
of their constituent solitons ]?]]. The purpose of our studies is to derive an explicit func- 
tional form for the next-to- leading-order (C(t -1 / 2 )) term of the effect of this interaction for 
the MNLSE. Although there have been several papers devoted to studying soliton solutions 
of the MNLSE ||, to the best of our knowledge, very little, if anything, is still known 
about its solution(s) for the class of non-reflectionless initial data. As is well known @, the 
investigation of multi-soliton solutions on the continuum background requires more sophis- 
ticated techniques, e.g., asymptotic, than, say, the algebraic methods for the construction 
of pure soliton solutions. An asymptotic investigation of the aforementioned solution can 
be divided into two stages: (i) the investigation of the continuum (solitonless) component 
of the solution |10|| ; and (ii) the inclusion of the soliton component via the application of 
a "dressing" procedure [|ll]] to the continuum background. The purpose of this paper is 
to carry out, systematically, stage (i) of the above-mentioned asymptotic paradigm (since 
this phase of the asymptotic procedure is rather technical and long in itself, the completed 
results for stage (ii) are the subject of a forthcoming article |l2|). The results obtained in 
this paper are formulated as Theorems 3.2 and 3.3. 

This paper is organised as follows. In Sec. 2, we transform (2) into a gauge-equivalent 
NLEE (see the first of (4)), and, using only some basic facts concerning its direct scattering 
theory |{13[|, pose the inverse scattering problem as a matrix oscillatory Riemann-Hilbert 
(RH) factorisation problem. In Sec. 3, we: (i) briefly review the Beals and Coifman for- 
mulation 1 14] for the solution of a RH problem on an oriented contour; (ii) give a general 
outline of the Deift and Zhou procedure j^] for obtaining the long-time asymptotics of 
the solution to the RH problem; and (iii) state our final results as Theorems 3.2 and 3.3. 
In Sec. 4, the initial oscillatory RH problem is formulated as an auxiliary RH problem on 
an augmented contour. In Sec. 5, the auxiliary RH problem is reformulated as an equiv- 
alent RH problem on a truncated contour. In Sec. 6, it is shown that, to leading order 
as i^+oo, the solution of the equivalent RH problem on the truncated contour is equal, 
modulo decreasing terms, to the solution of an explicitly solvable model RH problem on a 
contour which consists of the disjoint union of three crosses. Finally, in Sees. 7 and 8, the 
model RH problem is solved, and the asymptotic solution of related, auxiliary NLEEs and 
the MNLSE are obtained. 



2 The Riemann-Hilbert (RH) Problem 

We begin by introducing some notation which is used throughout the paper: (i) for D 
an unbounded domain of MuiM, let S(D) denote the Schwartz class on D, i.e., the class 
of smooth scalar- valued functions f(x) on D which together with all derivatives tend to 
zero faster than any positive power of | 1 for |x| — > oo; (ii) 03 = diag(l,— 1), a + = 

(° 1 

\0 

ro ad(<x3) T = m * aTm -* 3 . ( iv ) for k e { 1)2 } 5 £<£ x2 \d) = {F(X); A G D, ify(A) G C k {D), 
i,j G {1,2}}, where C k {D) = {/(A); A G D, \\f\\ cHD) = (JolfW^dXl) 1 ^ < 00}, and 

C& x2) (D) = {G(A); A G D, WGyWc-m) = sup|G 4J (A)| < oo,i,j G {1,2}}, with the 

xeD 









h- 


! 






and <Ti = er_ + 0"+; (iii) for a scalar w and a 2 x 2 matrix T, 
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norms taken as follows, ||(-)llr( 2x2 )mi = max ll(")ij'll£ n (D)) 71 G {1>2, oo}; (v) (•) denotes 

n \ ' *j's{i,2} 

complex conjugation of (•); and (vi) I denotes the 2x2 identity matrix. 

Proposition 2.1 ([]13[j) The necessary and sufficient condition for the compatibility of the 
following system of linear ODEs (the Lax pair) for arbitrary A G C, 

d x 9(x, t\ A) = U(x, t; A)\I>(x, t; A), d t ^>(x, t; A) = V(x, t; A)\&(x, t; A), (3) 

where 

U{x, t; A) = A(-iAcJ3 + Pa_ + Qa+) - %PQa 3 , 
V{x,t- 1 \) = 2\ 2 U{x,t-X)-i\{{d x P)a_-{d x Q)a + ) + ^{PQ) 2 ^(Pd x Q-Qd x P))a^ 

with P(x,t) = Q(x,t), is that Q(x,t) {resp. P(x,t)) satisfies the following NLEE, 

id t Q + d 2 Q + iQ 2 d x Q + \Q\Q\ A = (id t P - 8 2 X P + iP 2 d x P - ±P|P| 4 = 0). (4) 

Proof. (4) are the Probenius compatibility conditions for (3). ■ 

Proposition 2.2 Let Q(x,t) be a solution of the first of Eqs. (4). Then there exists a 
corresponding solution of system (3) such that *$>(x,t;0) is a diagonal matrix. 

Proof. For given Q(x,t), let ^f(x,t;X) be a solution of (3) which exists in accordance 
with Proposition 2.1. Setting A = in (3), we get that §(x, t; 0) = exp{-^ |Q(£ ; t)\ 2 d{,} 
■ IC, for some xq G R and non-degenerate matrix K, which is independent of x and t. The 
function ^(x, t; A) = ^(x, t; A)/C _1 is the solution of (3) which is diagonal at A = 0. ■ 

Proposition 2.3 ( [|l6f| ) Let Q(x,t) be a solution of the first of Eqs. (4) and ty(x,t; A) the 
corresponding solution of system (3) given in Proposition 2.2. Set ^f q (x, t; A) = ^ _1 (x, t; 0) 
•W(z,t;A). Then 

d x V q (x,t;X) = U q (x,t;\)V q {x,t;X), d t V q (x,t;X) = V q (x,t;X)^ q (x,t;X), (5) 

where 

U q (x,t; X) = -iX 2 a 3 + X(qa- + qa + ), (6) 
/ -2iX A -iX 2 qq 2X 3 q + iXd x q + Xqq 2 \ 

V q [X, t,X)-y 2A 3_ _ iXQx _ + x _ 2q 2 . a4 + . a2 _^ j , (7) 

g(s,t) = Q(x,t)((^- 1 (x,t;0)) 11 ) 2 , (8) 
is the "Kaup-Newell" Lax pair for the derivative non-linear Schrddinger equation, DNLSE 

id t q + d 2 x q-id x (\q\ 2 q)=0. (9) 

Proof. Differentiating ty q (x,t;X) = \& (a?, t; 0) i S(x, t; A) with respect to x and t and 
using the fact that 0) = exp{-^/* o i)| 2 a!£}, for some x G R, and ^(x,i;A) 

satisfy (3) for A = and A G C\{0}, respectively, defining q(x,t) as in (8), we get that 
^ q (x,t;X) satisfies (5), where U q (x, t; A) and V q (x,t;X) are given by (6) and (7): (9) is the 
compatibility condition for (5). ■ 
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Proposition 2.4 If q(x,t) is a solution of the DNLSE (Eq. (9)) such that q(x,0) G S(W), 
then u(x,t) = -^=exp{|(x - - x, |) satisfies the MNLSE (Eq. (2)) with u(x, 0) £ 

S(R). 

Proof. Direct substitution. ■ 

Definition 2.1 For Q(x,t) , as a function ofx, £ <S(M), define the vector functions ip^ (x,t; 
A) and ^(x, i; A) as i/ie Josi solutions of the first equation of system (3), namely, (d x — 
U (x, t; A))V> ± (x, t; A) = and (d x — U(x, t; X))cj) ± (x, t; A) = 0, with the following asymptotics, 

lim ^ + (s,t;A) = (0,l) T e ix2x+2iXit , lim ^~(x, t; A) = (1, ) T e - iA ^- 2iA4 *, 

x^+oo — x^+oo — 

lim + (x,t;A) = (l,0) T e- iA2a: - 2<A4 *, lim <jr (x, t; A) = (0, _i)T e iA 2 *+2iA 4 t j 

x— »— oo — x^ — 00 — 

where the superscripts ± mean $5(A 2 )^0, and T denotes transposition. 

Proposition 2.5 T/ie Josi solutions (Definition 2.1) have the following properties: (i) 

$^{x,t;\) = ^a ± (A,t)^(x,t;A) + 6^(A,t)0 ± (x,t;A), 3(A 2 ) = 0, 

wftere a ± (A,t) = VF^x, X),± ± (x, t; A)), ^(A,*) = T^(^(i, i; A),^ T (x, t; A)), and 
W(f,g) isthe Wronskian off andg, W(f,g) = hgi-fagi; (ii) if^(x,t;X) and c^(x,t;X) 
satisfy the second equation of system (3) as well, i.e., (dt — V(x,t; X))i/j^ z (x,t; X) = 
and (d t - V(x,t; X))^ ± (x,t; A) = 0, then a ± (X,t) = a ± (A,0) = a ± (X), and b ± (X,t) = 
b ± (X,0) e ±4ix4t = b ± {X )e ±AixH ; (in)forQ(X 2 ) = 0, a + (A)a~ (A)+6+(A)&- (A) = I, a+(A,t) = 
a-(X,t), b + (X,t) = -b-(X,t), a ± (-X) = a ± (X), b ± (X) = -6 ± (-A) ; and ^(P^) = 

-J? gj5Q g^T M fe with r±(X,t) = b±(X,t)/a±(X,t) r±{X)e±^ xH , and f = {A; 
Q=(A 2 ) = 0} (oriented as in Fig. 1); (iv) a ± (X) are analytically continuable to 9(A 2 )^0 with 





3(A) 


3(A 2 )<0 


'+ S(A 2 )>0 




+ 31(A) 


+ 

3(A 2 )>0 + ; 


- 3(A 2 )<0 



Figure 1: Continuous spectrum T. 

the following integral representation, ±ln(a ± (A)) = — jp ^^7~^~pj ~ J^/ (v,) /or 
A ->■ oo ; 9(A 2 ) > 0, a+(A, t) = 1 + 0(X~ 2 ), 

± + (x, t; \) e - ix2 *- 2ix4t = , l) T + 0(A" 2 ), 

<£+(x, t; X)e^ 2 ^ = (l, %^) T + 0(X~% 
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and, for A -> oo, 3(A 2 ) < 0, a" (A, t) = 1 + 0(A" 



Vf(x,£;A)e iA2a; + 2iA4 ' = (l,^) 

A)e- iA2 - 2iA4 * = - (%2 l) 1 + 0(A- 2 ) 



+ o(\- 

T 



Proof. See |§. ■ 

Remark 2.1. We now adopt a convention which is adhered to sensus strictu throughout 
the paper: for each segment of an oriented contour, according to the given orientation, the 
"+" side is to the left and the "— " side is to the right as one traverses the contour in the 
direction of the orientation; hence, (•)+ and (•)_ denote, respectively, the non-tangential 
limits (boundary values) of (•) on an oriented contour from the "+" (left) side and "— " 
(right) side (this should not be confused with the db superscripts). 

As a result of Propositions 2.1-2.4, for the asymptotic analysis of the solution to the 
MNLSE, it is enough to investigate asymptotically the solutions of (4) and the corresponding 
function ty(x,t;0) (Proposition 2.2): the main tool for this analysis is the following matrix 
RH factorisation problem. 

Lemma 2.1 Let Q(x,t), as a function of x , 6<S(R). Set ^>(x, t; A) = m(x, t; A) exp{— i(\ 2 x 
+ 2A 4 t)a 3 }. Define: (i) for 9(A 2 ) > 0, 



m(x, t; A) 



a+(X,t) %_KX,t,A) 



exp{i(A 2 x + 2A z 4)o 3 }; 



(ii) for 9(A 2 ) < 0, 



m(x, t; A) 



exp{i(A 2 x + 2A 4 t)a 3 }; 



and (in) r(A) = r (A). Then the 2x2 matrix-valued function m(x,t; A) (det(m(x, t; A)) = 1) 
solves the following RH problem: 

a. m(x, t; A) is holomorphic V A € C \ T; 

b. m(x, t; A) satisfies the following jump conditions, 

m + (x,t;X) = m-(x,t; A) exp{-i(A 2 £ + 2A 4 t)ad(cr 3 )}G(A), A G f, 

where 



G{\) 



1 -r(A)r(A) r(A) 
-r(I) 1 



r(A) G 5(f), and r(-A) = -r(A); 
c. as A — > oo, A G C\r, 

m(x,i;A) = I + C(A _1 ) 



Proof. Follows from the definitions of m(x,t;X) for 9(A ) < given in the Lemma and 
the properties of the Jost solutions and scattering data for (3) given in Definition 2.1 and 
Proposition 2.5: for details, see [19|. ■ 
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Lemma 2.2 Let IMI^oo^f) = su Pl r (^)l < 1- Then: (i) the RH problem formulated in 

Lemma 2.1 is uniquely solvable; (ii) ^(x,t;X) = m(x, t; A) exp{— i(X 2 x + 2A 4 t)o"3} is the 
solution of system (3) with 

Q(x, t) = 2i lim (Xm(x, t; A))i2 and P(x,t) = —2i lim (Xm(x,t; X)) 21] (10) 

A^oo A^oo 

(ra_j the functions Q(x,t) and P(x,t) defined by Eqs. (10) satisfy Eqs. (4), and 

q(x,t) = Q(x,t)((m- 1 (x,t;0)) 11 ) 2 (11) 

satisfies the DNLSE (Eq. (9)); and (iv) m(x,t;X) possesses the following symmetry reduc- 
tions, m(x, t; —A) = (73171(1, t; X)a^ and m(x, t; A) = a%m(x, t; A) 01 . 

Proof. The solvability of the RH problem Vi follows from the condition IM^oo^f) ^ 1 
and the discussions in (l^, [2^] (see, also, Zhou's skew Schwarz reflection invariant symmetry 
arguments [21, ^2|). The fact that q(x,t) defined by (11) satisfies the DNLSE follows from 
Proposition 2.3. ■ 

3 Solution Procedure and Summary of Results 

Prior to describing the Deift and Zhou procedure |l5| for obtaining the leading-order tem- 



poral asymptotics of (4), let us briefly review the Beals and Coifman formulation [14] for 
the solution of a RH problem on an oriented contour H. The RH problem (in the absence 
of a discrete spectrum) on H consists of finding a 2 x 2 matrix-valued function m(A) such 
that: (i) m(A) is holomorphic V A G C \ 3; (ii) m+(A) = m_(A)-u(A) VA G H, for some 
jump matrix u(A):H — » Mat(2,C); and (hi) as A — > 00, A G C \ E, m(A) — * I. Writing the 
jump matrix, v(X), in the following factorised form, v(X) = (I — w_(A)) _1 (I + w+(X)) V 
A G E, where w±(X) G C^ x2 \e) n £^ x2) (H), with \\w±\\ c (2x2)^ nC (2x2) ^ = ||w±|| £ (2x2) (H) 
+ Il u '±|| r (2x2),„ v are off-diagonal upper/lower triangular matrices, denote w = w + + w;_, 

and introduce the operator C w on C 2 (E) as follows, 

C w f = C + (fw-) + C-(fw + ), (12) 

where f££^ (E), and C±: £2 (S)— ^jC^ 2 (S) denote the Cauchy operators (bounded 
monomorphisms) , 



(C±/)(A) = lim jT 



A' e ± side of B 



/(0 # 

s (C - A') 2vri ' 



Theorem 3.1 (Q) J/ //(A) (= m+(A)(I - w+(A)) = m_(A)(I + w_(A))) G C {2x2) [E) + 
£oo X2 ^(3) solves the following linear singular integral equation, 

(Id - C w )n = I, (13) 

where Id is the identity operator on ' (E) + £^~> (E), then the solution of the RH 
problem for m(A) is 
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From Theorem 3.1 and (10), we obtain the following integral representation for P(x,t) 
(Q(x,t) = P(x,t)) in terms of the resolvent kernel of the linear singular integral equation 
(13), 

P(*,t) = -*([**, ^((Id-^.J^I)^)^^)]^) 
(see Sec. 5), which can be written in the following, simpler form, 

P(x,t) = - (Mx,t;i)) 22 {w^))2ie 2itm di : where 0(A) = 2A 2 (A 2 - 2A 2 ,), (15) 
7r Jr 

with real, first-order (distinct) stationary phase points, {0,±Ao} = {0, ±^^J— |} ({0,±Ao} 
= {A'; d A 0(A)| A=A , = 0}). 

At this point, the natural tendency would be to apply the method of steepest descents 



23] to the integral in (15); however, even though 0(A) is known and w-(\) can be char- 
acterised completely, no a priori explicit information regarding the analytical properties 
of the resolvent kernel, fi(x,t;X), is available. To eschew this manifest difficulty, Deift 



and Zhou [15|, based on the earlier, seminal work of Beals and Coifman [14|, introduced 
a novel non-linear analogue of the classical steepest-descent method and showed that, as 
t — > +oo, via a sequence of transformations which convert the original (full) RH problem 
(Lemma 2.1) into an equivalent, auxiliary RH problem with jump matrix v' of the form 
v' = fmodei + terror, where t'modei denotes the jump matrix for an explicitly solvable model 
RH problem and i> e rror contains only terms which are decreasing as t— >+oo, modulo some 
estimates, the solution of the original RH problem converges to the solution of the model 
RH problem. More precisely, the sequence of transformations which constitute the Deift 
and Zhou procedure [15| can be summarised as follows (even though the delineation of the 



procedure is general enough, the notation is specific to this paper): 



(i) based on the classical method of steepest descents |23j for decomposing the complex 
plane of the spectral parameter (A) according to the signature of K(z't0(A)), deform 
the original oscillatory RH problem to an auxiliary RH problem formulated on an 
(oriented) augmented contour X in such a way that the respective jump matrices 
on r C X and the finite "triangular" parts of X \ T, away from the neighbourhood 
of the real, first-order stationary phase points, {0, ±Ao}, converge (with respect to 
appropriately defined norms) to the identity (I), i.e., rewrite the original undulatory 
RH problem as an auxiliary RH problem on the paths of steepest descents with 
upper /lower triangular jump matrices whose structures depend on the signature of 
M(it0(\)) (see Sec. 4); 

(ii) reduce the contour X to a truncated contour X' which can be written as the disjoint 
union of three crosses, X' = X^/ U X#/ U X^/ (according to the number of real, first- 
order stationary phase points: see Sec. 5); 

(iii) analysing the higher order interaction between the three disjoint crosses, X^/, X#/, 
and Xc/, and noting that, as t — > ±oo, they are negligible (with respect to certain 
norms), separate out the contributions of the disjoint crosses in X' and show that the 
solution P(x,t) (Q(x,t) = P(x,t)) can be written as the linear superposition of the 
contributions of the various disjoint crosses (see Sec. 6); and 

(iv) localise the jump matrices of the most rapidly descented RH problems on the disjoint 
crosses to the neighbourhood of the real, first-order stationary phase points, {0, ±Aq}, 
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and, under suitable scalings of the spectral parameter, reduce the respective RH 
problems to RH problems on R with constant jump matrices which can be solved 
explicitly (see Sec. 7). 

In this paper, we prove that the application of the aforementioned procedure leads to 
the following results. 

Remark 3.1. To facilitate the reading of the results stated in Theorems 3.2 and 3.3 
(see below), as well as those which appear throughout the paper, the following preamble 
is necessary: (i) M G M>o is a fixed constant; (ii) the "symbols" ("notations") c and c 5 , 
respectively, which appear in the various error estimates are to be understood as follows, 
c = c(\o)eC co (R >M ), and c s = c s {\q) G S(R >m ); (hi) even though the "symbols" c and c s 
appearing in the various error estimates are not equal and should properly be denoted as 
Ci(Ao), £2(^0)5 etc., the simplified "notations" c and c s are retained throughout since the 
main concern in this paper is not their explicit (functional) Ao-dependence, but rather, the 
explicit class(es) to which they belong; and (iv) with the exception of Theorems 3.2 and 
3.3 (see below), the arguments of the "symbols" appearing in the error estimations will be 
suppressed. 

Remark 3.2. In Theorems 3.2 and 3.3 (see below), one should keep everywhere the 
upper signs as t^+co and the lower signs as t^—00. 

Theorem 3.2 For \ \r\ \ Coo ^ < 1, let m(x , t; X) be the solution of the RH problem formulated 
in Lemma 2.1 and Q(x,t) and q(x,t) be defined by Eqs. (10) and (11), respectively. Then 
as t^±oo and x^^foo such that Ao = \\J — f > M, 

Q(x, t) = exp{* (<S>± (Ac) + S ± (A ; t)) } + O , (16) 

q(x,t) = ^^exp{z (*±(Ao) + ^(A ;t)) } + o(&$±&), (17) 

where 

^) = -^ln(l-|r(z)| 2 ), (18) 

9+(z) = \ f z hV - z 2 | d\n(l - |r(^)| 2 ) - I / °° ln| M 2 + z 2 | dln(l + |rM| 2 ), (19) 

Q- (*) = I JT H^ 2 - A din(i - Km)I 2 ), (20) 

$ ± (A ;t) = 4A^^KAo)ln|t| ±argr(^(A ))+argr(Ao)T3KA )ln2, (21) 
$+(A ; t) = $+(A ; t) + f jf° Mlz^lQ^ _ 1 j- m(i + |rMl a ) d/x> (22) 

$ -(A ;t) = $-(*«,;*) + f J£ ^^'W , (23) 
and, as t—> ±00 and x—> ±00 such that A' = ^^/|>M, 

Q( x > *) = (^M) + @HK;t)) } + o(^^) , (24) 

= ^w$i^v{i (f (X' ) + &±(X> ; t)) } + O , (25) 

where 

v{iz) = -±\n(l + \r{iz)\ 2 ), (26) 

cf>+(z) = I J 2 ln|/i 2 - z 2 | dln(l + |r(^)| 2 ) - I / °° ln|^ 2 + z 2 | dln(l - |r(^)| 2 ), (27) 
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r(z) = \ JT ln|^ 2 - z 2 \ dln(l + |r(i/x)| 2 ), (28) 
e^A'o; t) = 4(A' ) 4 t t K**o) In |t| + e± ± arg r(*i/(»A{,)) + argr(iA' ) =F M^o) ln 2 , (29) 
8+(A' ; t) = §+(A' ; t) + £ / A ° MtB«l^ _ 1 iiM ^ ( 30 ) 
G -(A' ;t) = 0~(A o ;i) + $J% Mi+Wm*) dfij (31) 

e± = (2 ± 1)tt/4 ; and T(-) is i/ie gamma function 



Theorem 3.3 For | |r| | £OC( ,^ < 1 ; letm(x,t; A) be the solution of the RH problem formulated 
in Lemma 2.1 and u(x,t), the solution of the MNLSE, be given as in Proposition 2.4 in 
terms of the function q(x,t) given in Theorem 3.2. Then as A— > ±oo and x— >±oo such that 



*o=l/l(f-;)> M and f > |, sel >0 , 



= ^^W^) + (32) 



where 



A 2 , + — J t - KA ) In t - -j + arg r(ii/(A )) + arg r (A ) 

m 2 f^lnfl-Ku)! 2 ) , 2 A 00 ln(l + |r(uz)| 2 ) , , , 

- 2zy(A )ln2 + - / — i 1 Wl J dfi-- / — 1 v ^ n J dfj,, (33) 

7T Jo A* it Jo /X 

*~(Ao;t) = 2 (^o + ^) t + ^(A )ln|t| + ^ -argr(^(A )) + argr(A ) 

+ 2,(Ao)ln2 + ^r ln(1 " |rM|2) ^ (34) 

7T J\ Q fl 



and, as t— >±oo and x^=Foo or ±oo such that A = y ^(^ — f ) > M and § sGR>o, 



where 



e+(A' ;i) = 2^(A ) 2 -^yt-^A )lnt + ^ + argr(^(a )) + argr(fA ) 

/.f/M 2 /"^o ln(l + |r(iu)| 2 ) , 2 A 00 ln(l - |r(u)| 2 ) , , . 
- 2i/(iAo)ln2 + - / y 1 K J dfi- - — 1 Wl (36) 

7T Jo A* T JO A 4 



6-(Aj;t) = 2 ( (A'„) 2 - i- ) t+ f(iAj) ln |*| + ^ - aigr(«/(iAt)) + argr(iAi) 



+ 2 ,(,Ai)l„2 + ^r"' (1 + |r(i " )|2) ^. (37) 

7T JA' A* 

and !/(•), v(i-), and are given in Theorem 3.2, Eqs. (18)-(20) and (26)~(28). 
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Remark 3.3. The results presented in Theorem 3.3 give the asymptotic expansion of 
u(x,t) in the domain x/t ^ 1/s. In the case t — > ±00 and x/t = 1/s, u(x,t) can be written 
as follows, 

^,t) = ^ + o(^), 

where uq G C\{0} is some constant: the determination of uq requires special consideration 
and will be presented elsewhere. 

Remark 3.4- Some recent results concerning the asymptotic form of solutions as t 
— > ±00 of the DNLSE were written in terms of the so-called final states (the final-value 



problem) |25|]; however, unlike the present article, the authors in [25] did not treat the 
Cauchy initial-value problem. 

Remark 3.5. In the proof of Theorems 3.2 and 3.3, according to the above- formulated 
scheme, we omit the details of intermediate estimations which are analogous to those in 
15]: in corresponding places, we make exact references to the appropriate assertions in |i~5| ], 
and [19], where these details may be found; moreover, for Q(x,t), q(x,t), and u(x,t), one 
should consider actually four different cases, depending on the quadrant in the (x, i)-plane. 
In the subsequent sections, we present only the proof for the function Q(x,t) in the case 
t —* +00 and x —* — 00 such that Ao = \ \J—j > M: the results for the remaining domains 
of the (x, i)-plane are obtained in an analogous manner. The corresponding results for the 
functions q(x, t) and u(x, t) are obtained from Propositions 2.3 and 2.4. 



4 The Augmented RH Problem 

As explained in the previous section, we begin with the decomposition of the complex A- 
plane according to the signature of the real part of the phase of the conjugating exponential 
of the oscillatory RH problem (Lemma 2.1), 9t(it9(\)), where (9(A) = 2A 2 (A 2 - 2A§), A 2 , = 
-§ > 0, and t > (see Fig. 2). 





M(it0(A)) > // 


R(it0(A))>O 
» < 


1 < — '■> 




R(it0(A)) <0 


R(it0(A))<O 


-Ao _.--'6 


Ao] 


R(itfl(A)) >0 




K(it0(A)) < 0"'\ 



Figure 2: Signature graph of 3l(it0(\)). 



Remark 4-1- Note that, we have reoriented the contour T in accordance with the 
signature of §l(it9(\)) and the "sign" convention of Remark 2.1. 

The main purpose of this section is to reformulate the original RH problem (Lemma 
2.1) as an equivalent RH problem (see Lemma 4.1) on the augmented contour £ (see Fig. 3), 

S = LULuf, 
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where L = L° U L , and 



L° = {A;A = A + ^e^,ne (-00, ^2]} U {A; A = -A + ^e"T , u e (-00,^]}, 

L 1 = {A; A = ^fef,«G R}. 

Remark 4-2. Actually, for the following analysis, the augmented contour £ can be cho- 
sen in different ways, not necessarily consisting of straight lines: its important characteristic 
is the position of the contour L U L with respect to the lines 3t(itO(X)) = 0. 




Figure 3: Augmented contour S. 



In order to define the conjugation matrices on £ and exploit the analyses in [14, 
[l5[ |, we need to formulate two technical propositions: the first concerns the triangular 
factorisation of the conjugation matrices of the original RH problem (Lemma 2.1), and the 
second pertains to a special decomposition for the reflection coefficient, r(X). Write the 
jump matrices of the original RH problem (Lemma 2.1) in upper/lower triangular form for 
A G (—00, — Aq) U (+Aq, +00) as 



-ite(A)ad(cr 3 ) 



G(X) 



-M(X)ad(a 3 ) ( 1 r M 



1 







v 1 J\-r(X) 1 J 
and, in lower/upper triangular form for A G (— Aq,+Ao) U (— ioo,+ioo) as 



(38) 



G(A) 



-it0(A)ad(<73) 



1 



-r(A)(l -r(X)r(X)y 




1 



r(A)r(A)) 




;i-r(A)r(A)) 



-i 



1 r(A)(l -r(AMA))- 1 
1 



(39) 



To eliminate the diagonal matrix between the lower /upper triangular factors in (39), we, 
according to the scheme in |l5j, introduce the auxiliary function 5(X) which solves the 
following discontinuous scalar RH problem: 



S+(X) 



6-(X) = 5 (X), X G (-oo,-A ) U (+A ,+oo), 
<$_(A)(1 - r(X)r(X)), X G (-A , +A ) U (-too, +ioo), 
5(A) — > 1 as A — ► oo. 



(40) 
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Proposition 4.1 The unique solution of the above RH problem (40) can be written as 

5(A) = ^A^a)(A^a)) V+WeX-(A) e x+(A) e x-(A) } (41) 

where v = v{\§) is given in Theorem 3.2, Eq. (18), and 

X±{*)-2riJo ln {i-\r(\ )\*)~(jrXj> X±{X) - J± ioo (42) 

in Eq. (41), (A — n) " = exp{±i^ ln(A — //)}, where the principal branch of the logarithmic 
functions, ln(A — fi G {0, ±Ao}, /ias &een chosen, with the branch cuts along (— 00,/z), 
ln(A — yu) = ln|A — /i| + iarg(A — /x) ; arg(A — /i) G (— 7r, 7r); moreover, the function (5(A) 
possesses the following properties, 

5(A) = (stX))' 1 = (W^y 1 => 5(A) = <5(-A), (43) 
|5±(A)| 2 < (1 - sup^)r(/x)) ±1 < 00 V A G f, (44) 

MS? 

and IIW-))^!!^^) = SU P l(^(A)) ±1 | < 00. 

AeC 

Proof. It is well known that the scalar RH problem for <5(A) stated in (40) can be solved 
explicitly (see, for example, |26|, |2?]]): in our case, the solution is given in (41) and (42). 
As a consequence of (41) and (42), one proves the symmetry properties in (43): from (40), 
using (43), one deduces inequalities (44); therefore, from the maximum modulus principle 
and the fact that (5(A) has no zeros for A G C, one gets that (5(A) G £°°(C). ■ 

The conjugation matrices for the RH problem on the augmented contour £ should, 
of course, be written in terms of the matrix elements of the original RH problem on T 
(Lemma 2.1); but, since the reflection coefficient, r(A), does not, in general, have analytical 
continuation off T, we, following |l5| , decompose it as the sum of an analytically continuable 
part and a negligible non-analytic remainder. To formulate an exact result which we use 
later, let us define: (i) 



(x)= I r(A)(l-r(A)r(A))~ 1 , A G (-A , +A ) U (-ioo, +ioo), 
PK ' \ -r(A), AG(-oo,-Ao)U(+A ,+oo); 1 ; 

and (ii) the contour = Zii U Li, where 

K ' 5 5 <5' 

L^={\;\ = \ Q + ^e 3 -^,u£ (5, v^]} U {A; A = -A + e"T , u G (?,V5]}, 
Li={A;A = ^feX,nG [-y/2, -5) U (6, \/2]}. 

Proposition 4.2 For each I G Z>i, there exists a decomposition of the function p(X) in 
Eq. (45), 

p(X) = hi{X) + R(\) + hu(X), AGf, (46) 
such that hi(X) is analytic on T (generally, it has no analytic continuation offT), R(X) is 
a piecewise-rational function such that ^f^P = d3 j\P > — 3 — + 1, where 13 



xeu 



= {0, ±Ao}, and hu(X) has an analytic continuation to L; moreover, in the domain X$>M, 
the following estimates are valid as t— >+oo, 

|e-*" (A) WI < CTW ' AGf ' l^ 2 ^ (A) ^(A)|< (1+|A| f )(A ^ , AGL, 
\e- 2it6 ^ R(X)\ < cexp{-2)ffit}, X G Lj, 
and 5 G M>o is sufficiently small. 
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Proof. Rather technical: proceed analogously as in the proof of Proposition 1.92 in 
|j~5| ] by expanding p(X) in terms of a rational polynomial approximation in the neigh- 
bourhood of the real, first-order stationary phase points, {0, ±Ao}, and show that, for 
I G Z>!, \e- 2it9 Whj(X)\ < (i^W , A G f, \e- me ^h H {\)\ < A G L, and 

\e~ 2ite ^ R(\)\ < cexp{-iAo(?|x|}, A G L^. Using the relation \x\ = 4Agf, one obtains 
the result stated in the Proposition: for details, see [19|. ■ 

Lemma 4.1 Let m(x,t;X) be the solution of the RH problem formulated in Lemma 2.1. 
Set m A (x,t; X) = m(x,t; A)(A(A)) _1 ; where 

A(A) = (5(A))"», (47) 

and (5(A) is given in Proposition 4-1- Define 

m A (x,t;X), A G Q x U fl 2 U fl 3 U 4j 
m^ac,*; A) = <( m A (x, t; A)(I - (w^)^)" 1 , A G 5 U 6 U 7 U 8 U 9 U fi l0 , 

m A (x, t; A) (I + (w+jx.t^r 1 , A G On U O i2 U Oi 3 U 14 U O i5 U Oi 6 . 

Then (x , t; X) solves the following (augmented) RH problem on S, 

mJ_(z,t;A) =ml(a;,t;A)4, t ,«(A), A G E, ^ 
m"(ic, i; A) — ► I as A — > oo, 

where 

( (l-(w^) xt5 y 1 (l+(w + ) xt5 ), AGf, 
<^(A) = (I-(^i)^)- 1 (I+(4)^)= AGL,_ (49) 

I (I-(^)xA*)- 1 , AGL, 

and 

(w^' a) ) x ,t,s = (5±(X)r d ^exp{-ite(X) a d(a 3 )}w^ a \ (50) 
< = hi(X)a+, wl = {h H (X) + R(X))a+, (51) 



w° = -/u(A)o-_, uj« =-(^j/(A) + fl(A))<7_. (52) 

Proof. In terms of the function m A (x,t; A) defined in the Lemma, the original oscilla- 
tory RH problem (Lemma 2.1) can be rewritten in the following form, 

m A (x,t;X) = m A (x,t; A) (I - (w-) X}t ,s)~ 1 0- + (w + ) X}t ,s), A G T, ^ 
m A (x, t; A) — > I as A — > oo, 

where 

(u>±)*,M = (<5±(A)) ad ^)exp{-^(A)ad(a3)}y; ± , (54) 

and 



tu + = p(A)<7 + , = — p(A)cr_. (55) 

Defining m"(x, i; A) as in the Lemma, one arrives, as a consequence of Proposition 4.2, (47), 
and (50)-(55), to the RH problem stated in the Lemma. ■ 
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5 RH Problem on the Truncated Contour 



The goal of this section is to get rid (asymptotically) of the contribution of the func- 
tions hj(X) and hjj(\) (Proposition 4.2) to the conjugation matrices of the augmented RH 
problem: simultaneously, it is proved that, as t — > +00, under the condition Ao > M, the 
contribution to the asymptotics of the function P(x,t) (Q(x,t) =P(x,t)) coming from T 
and the finite "triangular" segments L^-and L~, respectively, are negligible; hence, we derive 
the RH problem on the truncated contour £' (see Fig. 4), 




Figure 4: Truncated contour S'. 



E' = E\(i ? UL ? uf). 
Instead of (10), from Lemma 4.1, one can write 

P(x,t) =nim(A[fT 3 ,m a (x,t;A)]) 2 i, Q(x,t) = P(x,t); (56) 

whence, following the Beals and Coifman formulation (Theorem 3.1), by setting 

w l,t,5 = ( w +)x,t,S + («>- )*,*,*• 

where {w^.) x ,t t S are given in (49)— (52), one gets the following integral representation for 
P(x,t), 

p( x ,t) = -i Qf to, ((M - e^J^m^tM^i)^ ' (57) 

where (Id — C a ) is invertible as an operator in C^^i^) + ^c» X2 ' ) (S) Q. To reduce 

x ,t,S 

the RH problem on S to the one on the contour we need to estimate corrections which 
arise as a result of the corresponding reduction of the integration contour in (57). To do 
this, we decompose ur x t s as 

w l,t,6 =w e + w', w e = w a +w b + w c : (58) 

(i) w a = w x 1 5 |p is the restriction of vf x t $ to T and consists of the contribution to w* Bt& 
from hj(X) and hj(\); (ii) w b has support onlUL and consists of the contribution to w x t s 
from /177(A) and hu(X); and (hi) w c has support on L-^UL^ and consists of the contribution 
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to w xtS from R(X) and R(X). The main idea of the estimation of the integral in (57) will 
be (based on Proposition 4.2) to show that, as t — > +00 (x/t ~ 0(1)), it/| s \jy — ► and 
w e — >0, in the sense of appropriately defined operator norms, and, therefore, to "lump" the 
contribution to %ir x t s from R(X) and R(X) into the factor w', which is supported on and 
show that it encapsulates the leading order asymptotics. 

Lemma 5.1 For arbitrary I € Z>i and sufficiently small 5 E M>o ; as £ — > +00 sttc/i t/iat 
A >M, 

ll/m^ll << — I Lm&I I ^ — I L11C 1 1 ^ 2An<5^£ 

II™ H4 2x2) (f)^(W l |w H4 2x2) M-(W 11^ Il4 2x2) (w-- ' 

||^ , || £ (2X2) (S) <^| T7 4, H W, H £ (2x2) (s) <^|p, 

where k G {1, 2, 00}. 

Proof. Consequence of Proposition 4.2, the following estimates, 

(«4)., M , 4 je £[ M) (E)ne 2) (s), fee {1,2}, 

and analogous calculations as in Lemma 2.13 of [15|: for details, see jnj. ■ 
Definition 5.1 Denote by Jx(-) the space of bounded linear operators acting in C 2 (•)• 



Lemma 5.2 ,4s i -> +00 suc/i i/taf A > M, (Id - C a ) _1 G Af(E) 44> (Id - C^y 1 € 
A/-(E). 



Proof. Consequence of the following inequality, IIC j — Cu/|U/7S) < c||w e || (2x2)™ , 
the fact that ||w e || (2x2), ,(< c\x\^ 1 ) <c(Xnt)~ l <c ((58) and Lemma 5.1), and the second 

Loo (2^) 

resolvent identity. ■ 



Remark 5.1. Actually, the operator (Id — C w >) 1 acts in /^^(^O + £oo X2 ^(E), and 
the operator C j acts from £[> 2x2 ^(£) + £oo X2 ^(E) into £*2 X2 \t,): we consider their 

x,t ,5 

restrictions to £^ 2x2 ^(E). 

Proposition 5.1 If (Id — C w >)~ 1 eAA(E), then for arbitrary ZeZ>i, as i— >+oo smc/i i/iai 
A >M, 

P(M) _ -,: (/ E fe((!d-0)-'.)(5K(0]i) !i+ o( i | y ). (59) 

Proof. From the second resolvent identity, one can derive the following expression (see 
Eq. (2.27) in [H 



hm-cj )- l i)(o<*. 5 (oll=/ s ((id-^)- 1 i)(e)^(e)^+/+^+m+^ 

where 
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The terms /, II, III, and IV have, respectively, for I G Z>i, the following estimate, 
0(c\x\~ l ) (hence, using the relation |x|=4A^, C(c(A^) _i )). Let us prove, for example, the 
last estimate: 

llVl^Uld-C^C^ild-C^ )~\C wt 1)1^(2X2)^11^^(2X2)^ 

x,t,S x,t,S 2 V / 2 V / 

1 1 tt 

x,t,S x,t,S ^2 \ J ^2 \ / 

< C| I C TO e| | | C7j ^ ^ 1 1 I £ (2 x 2) (s) | | Wj )tj a | | £ (2 x 2) (s) < c||u; e || £ (2x2) (s) (||^ ti( 5|| £ (2x2) {E) ) 2 . 

From Lemma 5.1, the arithmetic mean inequality, and the Cauchy-Schwarz inequality, one 
shows that, 

(K^lzf-W 2 - 4 (ll wa ll4 2 x 2 >(r)) 2 + 4 (ll^ll4 2 -)(LuI)) 2 

+ 4(|K||^ ( ^) a + 4(||«/|| jC?xa)(I g) 2 <c; 

whence, estimate \IV\ follows from the fact that | \w e \ L(2x2) <c(Ao£) _/ ((58) and Lemma 
5.1). ■ 

Let us now show that, in the sense of appropriately defined operator norms, one may 
always choose to delete (or add) a portion of a contour(s) on which the jump is l+0( n § .A , 

ZsZ>i and arbitrarily large, and without altering the RH problem in the operator sense. 

Definition 5.2 (p5|) Let: (%) R^> : £^ 2x 2 -*(£) — > C^ 2x2 \y,') denote the restriction map; 
(ii) I E ^ E :4 2X2) (S') -» 4 2x2) (£) ^noie tfie embedding; (Hi) C%:£%* 2 \n) -» 4 2><2) ( s ) 
denote the operator in Eq. (12) with w *-* w' ; (iv) C^':4 2X2) (S') -> 4 2X2) ( S ') de ^ e 
operator in Eq. (12) with w <-> u/|s'; (wj C^,: j0^ x2 \~E') — > 4 2><2 ' ) (^') denote the restriction 
of C^, to £2 v^')> and (vi) Id S / and Id s denote, respectively, the identity operators on 
4 2X2) (S') anrf4 2X2) (S). 



Lemma 5.3 



(Ms' - O -1 = ^'(Ids " C^r'ls^E, (60) 
(Ms " C*)- 1 = Ids + C5(Ids' " Cg)" 1 ^/. 

Proof. See Lemma 2.56 in [jig] . ■ 

Proposition 5.2 // (Id s — C^ ; ) € JV(S), i/ien /or arbitrary ZeZ>i, as t^+cx) sue/t i/iai 
A >M, 

P(M) _ -,: (£ K (dfc " cS)-W«))^) a + «>( I ^ r ). (6!) 

Proof. The boundedness of 1 1 (Idy; — C^')" 1 ! Ia^(S') follows from the assertion of the 
Lemma and identity (60): the remainder is a consequence of Proposition 5.1. ■ 

From Proposition 5.2, it is clear that the asymptotic expansion (to 0{ f ), ZeZ>i) 

can be constructed by means of the following RH problem on the contour 

(62) 



m%(x,t;\) = m?{x,t;\)v%' tS (\), A G 



■mP' (x, t; A) — ► I as A — > oo, A 6 C \ S', 
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wh ere v^ s (X) = (I - (w?)^)- 1 ^ + (w+)z,t,s), with 

«)*,*,* = (<5(A)) ad (-3) e -^(A)ad(a3) jR(A)(7+) (^^ = 0, A G L\Lg, (63) 
(«>?')*,*,* = 0, (^%, t , 5 = -( ( 5(A)) ad ^) e -^(A)ad(.3)7p) (J _ ) Ael\^. (64) 

Denote {vP')^ = (to+)s,t,<s + (u>E') x ,t,<s, so that (w E ' ) X)t)( j = u/|s'; then, according to 
Theorem 3.1, the solution of RH problem (62) has the following integral representation, 

m (x,t;A)=I + y E# ^ — , AEC\E, 

where ^'(x,t;X) = (Id s , - Cg)" 1 ! 

Remark 5.2. In (64), i?(A) is the same piecewise-rational function R(X) appearing in 
(63), except with the complex conjugated coefficients. 

6 RH Problem on the Disjoint Crosses 

In this section, we make a further simplification of the RH problem on the truncated contour 
£' by reducing it to the one which is stated on the three disjoint crosses, S^/, £#/, and T,q>, 
and prove that the leading term of the asymptotic expansion for P(x, t) (Proposition 5.2, 
(61)) can be written as the sum of three terms corresponding to the solutions of three 
auxiliary RH problems, each of which is set on one of the crosses; moreover, the solution of 
the latter RH problem can be presented in terms of an exactly solvable model matrix RH 
problem, which is studied in the next section. At the end of this section, we also prove the 
basic bound on (Ids' — C^') 1 (Proposition 5.2). 

Let us prepare the notations which are needed for exact formulations. Write £' as the 
disjoint union of the three crosses, £4/, £b' ; an d £c"i extend the contours £4/, £#/, and 
£(7/ (with orientations unchanged) to the following ones, 

t A , = {A; A = -A + Y e±f . u G M i> £ b' = {A; A = A + ^e ±2 f\ u G R}, 

g c , = {a ; A=^ ±? ,«el}, 

and define by £4, £5, and £c, respectively, the contours {A; A = ^^e^^ , u G M} oriented 
inward as in T,a> and £4', inward as in £b' and £b', and inward/outward as in T,c and 
£c- 

For /c G {^4, B, C}, introduce the following operators, 

N k :tf(%)^£\^ k ), /(A)~(tf fc /)(A) = /(A fc + £fc ), (65) 

where 

A A = -A , Xb = Xq, A c = 0, e A = e B = XiWXlty 1 / 2 , e c = A(8A^)-V2. ( 66 ) 

Considering the action of the operators iV^ on 5(X)e~ tt9 ( x \ we find that, for k G {A 5, C}, 
I A = (-00, Aa), /b = (As, +00), and I c = (-A , +A ), 

JV fc {5(A)e- i «W} = ^(A), K(A)G4, 

where, as a result of the second of (15), (41), and (42), 

<5? = (16A#)-$exp(2iAgt+ £ Xm(W , *€{A,B}, $ = expj £ Xro (0)l, (67) 

!. m£M I ImeA^ I 
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W = (X Xl%P yV exp (-T (1 + + ^,OCm(Ai + ei) - Xm(Ai))|, (68) 
^(A) = (e^-Agrexp{f (1 - + E (Xm(ec) - Xm(0))} , (69) 
with X = {A, B, +, -}, and 

Xk(Xi) = ^J Xk ^HM*)l^k-)> X ± (AO = ^ooH^-A^ln(l + |rM| 2 ), (70) 
Xfc(O) = h Io k HmI dln(l " l^)| 2 ), X±(0) = j£x, 

Set 

A2 = W», (7i) 

and let A° denote right multiplication by A£: 

A°*=*A°. (72) 



Denote 



^ | f >^ ^f' md = { of 'a 6 a la (73) 

According to this, 

(tO.A*= E = E C£,= E C^,': (74) 

fee{A,B,C} 1 ;a! '*' 4 fe6{A,B,C} fcG{A,B,C} 

hereafter, we do not introduce a special notation for w k \v ■ Let us prove some technical 
results concerning the operators C E K and c5K . 
Proposition 6.1 For k G {A, B, C} ; 

C*|' = (^)- 1 (A°)- 1 c2(A0)iV fe , w k ee (Al)-HN k w*')(A°), (75) 

u>/iere 



C^I £ ( 2 x 2 ) ( ^ ) = -G+(-(^ + A fe )(^(A))- 2 f 7_)), 



(76) 



i/ere, i/ie rays are defined as follows, 

U = {A; A = ^(lGAg^^e-x, u G (_ £) +0 o)}, Z G {A, B}, 
L c = {A; A = *f (8A 2 t)V2 e f j u G R} ? 

so i/iat Sfc/ = Lfc U Lfc. 

Proof. We consider the case k = B: the cases k = A and C follow in an analo- 
gous manner. Since, from the first of (67), = 1, it follows from the definition of 
the operator A° B in (72) that (A'g)t = (A^) -1 , where f denotes Hermitian conjugation 
(A° B is a unitary operator). Changing variables ((65)), recalling (12), using (71) and 
(72), as well as the unitarity of A° B , one obtains (75), where C^§ = C^ A o )-i^n b w b ')(A ) 
= C + (iA° B )- 1 (N B iu B ')(A B ))+C-(-(A B )- 1 (N B wf ){A%)). Using (63) and (64), one shows 
that ((A° J )- 1 (iVB^')(A0 J ))(A) = on and {(A° B )- 1 (N B wE')(A° B ))(X) = on L B , so 
that (76) are valid. ■ 
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Lemma 6.1 Let k G (0,1), p(B) = -p(A) = I, p(C) = -Ag/(2A), sgn(A) = sgn(5) = 
— sgn(C) = l. T/ien V AG Lfc C Sfc; as t ^ +cc such that Xo > M , 

\R(e k + A fc )(^(A))- 2 -i?(A±)(2Ap(fc))- 2 - e is s n ( fe ) A2 | < £©^-2^*^ (??) 



and VAeLfcCSfc, 



Iflfofc + A fc )(^.(A)) 2 - i?(A±)(2Ap(fe)) 2 -e- iS s n ( fc ) A2 | < £(^L* e -2^ 2 , 



(78) 



lim J R(A)=r(A )(l-|r(A )| 2 )- 1 ; 

;a(A) f Aq 



where (resp. L&), A; G {A, 7?,C}, are defined in Proposition 6.1, u 6 (— e, +oo) ; wi/t 
< e < v/2, i?(A+) = lim i?(A) = -r(A ), i?(A„ 

»i(A)J.Ao 

i?(-A+) = lim R(X) = -r(Ao)(l - KAo)! 2 )" 1 , «("A ) = Km fl(A) = r(A ), 

iK^AJJ.— Ao iK^AJT — Ao 

i?(0 ± ) = since r(0) = r(iO)=0, c(B) = <? , c(A) = c' s , andc(C) = c. 

Proof. We prove inequality (77) for k = B: the other results follow in an analogous 
manner. By using (68), for A G Lb C S_b and k G (0, 1), write 



i?(A + e B )(o B (A))- 2 - i?(A±)(2A)- 2w e^ = e^(7 + 77 + 777), 



2ivi\ 2 



where 



(2A)- 2w e^ 1 -f »\R(\o + e B ) - R(>£)), 



II = (2X)- 2iu e< 1 -^ x2 R(X + e B 
III = (2A)- 2 - e < 1 -f) A2 J R(Ao + e B 

and, in the last formula, Z = zA 2 f |^ + ^ 

| exp{^}| = exp{-2 K A(^ 2 }, v 
and 777 can be estimated in the following way. 

|7| < \{2Xy 2iv 



Ao(A + £ B /2) 
(A + e B ) 2 
A (Ao + £ B /2; 
{Xq + eb) 2 



-2iv 



- 1 



-2iv 



2 £ (Xm(A + £B) -Xm(Ao)). Note that 
exp{^y-}| = exp{— 2KAgtti 2 }, which gives the exponential factor in (77). The terms 7, 77, 



= i(l-f)A 2 



\eb\ sup 
Aeli' 



9 A «(A) 



< e 2 



^lAle-^-DW 2 



16Agt 



£°°(Lb) 



< 



where c is independent of A. 



1771 < e 2 e" 



-(i-f)|A| 2 



i? 



where 



77i 



/ Ao(A u + £B /2) y 
V (Ao+es) 2 J 



-2ii/ 



f A„ \ 
\\ Q +e B /2) 



-2iv 



Ih 



\Xo+£ B J 



-Aiv 



To estimate |772|, one proceeds as follows: 

fr (l- f )|Ah // , <e -(i_ S )| A |2 



< 4^( 1 -f)l A l 2 sup{|^- 1 |;e = 1 + G [0,l]}J^i < ^t" 1 / 2 , 
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since \^ 2%u l \ < v2 exp{— 2i> arg(£)} and — ^ < arg(£) < f: in this estimation, one uses 
the fact that |£| > l/\/2 and < v < y max = — ^ ln(l — supr(/Z)r(//)) < oo. Since the first 

Me? 

term on the right-hand side of the equation for II\ is bounded, namely, 



A (Aq + £b/2) 

(A + e B Y 
one gets an analogous estimate for |i7i|. 



-2iu 



< e 2 



II/II < e" 



< ce 



R 



e -(l-f)|A| 2 gup 



C°°(L B ) 
~ ^75 )|A| \Z\ < ce 



0<s<l 



■e s * 



f-i73)l A l 2 f 



< ce -(i-f)|A| 2 | Z | ™- 1^1 



sup |e 

0<s<l 



A4 + o 



c 5 lni 



16Afc 



where one uses the boundedness of exp{ E (xm(Ao + sb) — Xm(Ao))}, which follows from 

(44), (67), and (68). The term 0(c s {\^t)~ 1 / 2 hit) appears in |///| due to the estimation 
Xz(Ao + £b) — X/(Ao)> ^ £ {A,B}, which can be obtained by using the Lipschitz property 



of the function In 



i-Hm)I 2 

l-KAo)|* 



I, \fi\ <Ao, and integrating by parts the first of (70): analogous 
estimations for x±(Ao + £b) — X±(Ao) are 0(c(Agi) _1 / 2 ). ■ 



Proposition 6.2 For general operators C^' k , , k £ {1,2, .. . ,N}, if (Id S / — C s ' k ,) 1 exist, 



then 



(id S '+ e ^(id^-c^:,)- 1 )^,- e cj;,) 



Ka<N 



l<f3<N 

Ms' - E E (i-M(Me'-C')~ 1c7 S' c S" 

l<a<Af 1</3<AT 



and 



(Id E ,- E cS)(M S '+ E c^dd^-cS)- 1 ) 



1</3<N 



Ka<N 



Ms' - E E {i-5 aP )c* a ,c^,{Uv-c^,)-\ 

l<a<Nl<(3<N 

where 5 a p is the Kronecker delta. 

Proof. Assumption of existence of general operators (Id S / — C^ fc/ ) _1 , k £ {1,2, ... ,N}, 
induction, and a straightforward application of the second resolvent identity. ■ 



Lemma 6.2 For a / j3 G {A',B',C'}, as t^+oo such that Xq>M, 



\C, 



l/^S' ^E' 



< 



Proof. Analogous to Lemma 3.5 in [15]. 



Lemma 6.3 //, for k G {A,B,C}, (Id E , - C S £,') _1 G AA(S fc /), t/ien as t^+cx> stic/i that 
X >M, 

p(*,t) = -i e (/ EA j^,((M Sfe ,-c'5;)- i i)(e)^'(e)]S) 21 + (^)- ( 79 ) 

k€{A,B,C} K k W /21 
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Proof. Using Proposition 6.2 and the second resolvent identity, one writes 

(Ms'- E CjlO-^D^+D^W^-Es,)- 1 ^,, (80 ) 

ke{A,B,C} 



where 



Djy = Ids' + £ C*V (Ids/ " C\, ) " 1 , 

(Ids' - ^ 



Es' - E (i - W^^dfc - c^) 

a,/3£{A,B,C} 



Since, from Lemma 5.1, ||C , ^,||^(s') < Il ti;fc 'll£(2x2)^ s ^ < ^2^171 < c, and, by assumption 
1 1 (Ms' ~~ C^1') _1 |Ia/'(S') < °°> v i a Lemma 6.2, one finds that, as such that Ao>M, 

HI>s'lk(S') < c, and ||(Id s , - E S / ) _1 1 ^(^z) < c. 

Taking account of the second resolvent identity, one also finds that, 

||E S J|| ( 2x2), SM < E (1 - ^/3)l|C^ Q 'C^/3' I ll/ : ( 2x2 )(-S'-l 

+ e (1 - ^)iicj:,cj;,m s0 n(id s , - c^iw^c^Mi^ 

a,l3€{A,B,C} ww w w L 2 ^ J 



Noting from Lemma 5.1 that, 



l|C,,,fe'I|l r (2x2),„ M < \\W || (2x2),™, < 



one finds, using Lemma 6.2, and applying once more the second resolvent identity, 

||E S ,I|| £ (2x2) (s/) < 3^74. 

Now, from the Cauchy-Schwarz inequality and Lemma 5.1, 

||E E ,U> S || £ (2X2) (S;) < ||E S ,I|| £ ( 2X 2) (E , ) ||W S II £ (2X2) (S , ) < Ao |/4 < ^ ! 



hence, recalling (80), (Ids' ~ E C S 1') _1 £ AA(S'). From Proposition 5.2 and the 

k£{A,B,C} W 
above estimates as such that Xq>M, 

J s ,((Id s , - ^')- 1 I)(0^ E '(0^ = J S ,(D S ,L> S '(£R + O(^) + 0(jfa), (81) 

for arbitrary / G Z>i. Recalling that w s ' = X) the integral on the right-hand side 

ke{A,B,C} 

of (81) can be written as follows: 

/ (D s ,iv s '(eR= //Ms' E E ^((ids'-^l')- 1 !)^')^)^. 

S S ' ke{A,B,C} a,(3£{A,B,C} 

(82) 

To estimate the right-hand side of (82), consider, say, the following integral, J S /(C S A , 
•(Id^-C^)- 1 !)^'^: 

i/ s ,(cj;,(id s , -c^r^io^mi < i/s'(^i)(o^ B '(o^i 
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MiAc^m^-c^c^m^'imi <\iu 

((Id s ,-C S ' )- 1 C S ' I)(|,)W A '(D) , D , (2X2) ll» S 'll (2X2) 

+1 J U tFS UJ^I < 2^dist(s A ,,E B ,) 

ii^-gS)" lc S I ''4»><»)(^J'^ , ''4»><»)( SA< )''"' fl, '' £ ('^)( Sfl< ) 

' 2tt dist(S A ,,S s ,) 

< 1 1 \W A ' || £ (2X2) I | £ (2x2) (Eb;) 

+^ 1 1 (Ids' - C^l' )~ X I Ia^(s^) I IC'^' 1 ! Ufx 2 ) 1 1 I x2) (s^) " 1 Irf x2) ' 

since 2tt dist(S J 4/, Eg/) > Aq/c. Hence, by Lemma 6.2 and the assumption that ||(Id S / — 
C^l') _1 |lAf(£') < oOj one gets, for a ^ j3 £ {A, B,C}, the following estimation: 

therefore, applying the second resolvent identity to the right-hand side of (82), one finds 
that, 

S ' k&{A,B,C} S ' 

for arbitrary / gZ>i. Now, from Definition 5.2 and Lemma 5.3 (identity (60)), 

(Id Sfc , - C^T 1 = i? Sfc ,(Id E , - C^)" 1 ^^; (84) 
hence, substituting identity (84) into (83), and recalling (80) and (81), the proof is complete. 



+ Ol S V (83) 



Lemma 6.4 For k e {A, B, C}, (Id S , - C^v) -1 G -^( s 



Remark 6.1. This result was proved in [15]: here, we briefly reproduce this proof since 



we need a model RH problem which arises in it in order to obtain the explicit asymptotic 
formulae presented in Theorem 3.2. 

Proof. Consider, say, the case k = B: the cases k = A and C follow in an analogous 
manner. Define the function (see Fig. 5(a)) 

luf (A) = ({A%)-HN B w¥)(A%))(\) = ±i? ± (A )(2A) ±2 -e^ 2 C T ± , 

where w± (\) is given in (73), and 



Now, defining as usual w B ° = w+° + w B ° , and using Lemma 6.1, one finds that, 

\\w B - u^Ha = ||(A^)- 1 (JVBi5 B, )(A^) - (A b )-^ b zI; b0 )(A° )|| N < 4§* (85) 
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V 4 



V 2 
B 




E s (b) (c) 

Figure 5: (a) T, B ; (b) and (c) S e:r = S B , r U M. 



(2 x 2) 

where H = £p (£b), pG {1, 2, oo}. Hence, as i— > +oo such that Aq > M, 



\\ C J ~ C w Bo\\c<?* a \v B ) ~ 



(86) 



and, consequently (from the second resolvent identity), one sees that, for sufficiently large 
t, (Id SB - Cl B B0 y l G AA(S B ) (Id Ss - C^iY 1 G Af(T, B ). Reorient S B as in Fig. 5(b) 

and define w B,r (X) = w+' r (\) + w B ' r (X), where 



u>£' r (A) = w B "(X), 5ft(A)>0, and w± r (X) = -w B "(X), »(A) < 



so that one can consider the operator C^f;^ = C + (-w B ' r ) + C_(-w+' f ), where the Cauchy 
operators are now taken with respect to Extend T,B, r — ► ^ex = £,B,r u R, with the 

orientation in Fig. 5(c), and set 



w ex (X) = 



w B ' r (X), AG£ B , r , 

0, A G £ ea; \ S_B jr . 



From Lemma 5.3, it follows that (Id_s ea; - C^) -1 G A/"(S e 



(Ids 



Sfl.r-, — 1 



€ 



Af(£ Bjr ) (Ids B - C^o)" 1 G AA(S B ). Define a piecewise analytic 2x2 matrix-valued 
function </>(A) as follows, 

4>(X) = { 



(2A)- < "°3, AgO|UO§, 
^-^(I-tuf^A))- 1 , A€fifUJ21, 
(2A) -iI " 73 (I + to^f^A)) -1 , Aefi§ung, 



(87) 



where I± tu^A) are, respectively, the analytic continuation of I± w ex (X) into J7f, 
and 0|. From the estimates above for w±°(X), one finds that, 

wff (A) — > as A — ► oo, A G Of U Q|, 
i^(A)^0 as A^oo, AGOfung. 

For A G S e:c , set V e ^{X) = 0_(A)(I - u/I*)- 1 ^ + u»f )0+ 1 (A) = (I - wff'*)- 1 ^ + 
where (f>±(X) are the non-tangential limits of <p(X) as A approaches T, ex from the "±" side, 
respectively (Remark 2.1). For A G E B , r , V e ^{X) = I, and, for A G T, ex \ E B , r , V e '^{X) = 
(f>-(X)(f>'T 1 (X), whence, one gets that, 



V e ^{X) = exp {-^ad(a 3 )} v(X ), 5ft(A) < 0, 
V e ^{X) = exp {-if ad(a 3 )} v(A ), »(A) > 0, 



(88) 
(89) 
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where 



«(A„) = ( 1 -Jgf det W A )) = l. (90) 



Hence, the jump matrix, V e, ^{\), can be characterised as follows: 
V e ^{\) = 



I, A G S^r, 

exp{-^ad(cr 3 )}w(A ), A G S e:E \ E B , r - 



On IR, one has that, 

V e ^{\) = (I - u;^' )" 1 ^ + w^*) = (I - f(V)e a V)(I + r(A )e- u2 c7 + ). 

Set C ey(j) = C + {-w e * ,4 ') + C-(-w+' <t> ) as the associated operator on S e:c , with w ex ^ = w+ ^ + 

w e*,<P_ By Lemma 5.3, the boundedness of C e ^, i.e., ||(Ids ex ~ C^S.^)" 1 ! k(E e ,), follows 

from the boundedness of the operator C w ex,<t>\ : C^* 2 ^ (IR) — ► C^* 2 ^ (IR) associated with the 

IR 

restriction of u; 6 ^ to IR; but HC^e*,^ ^II^R) < sup |r(A ) exp{-iA 2 }| < ||r|| £oo( ^ < 1, 
and hence, by the second resolvent identity, 

11(1% - C^i^)- 1 !!^ < (1 - llrll^^)" 1 < oo 
(Id M - C^^)" 1 G AA(M) =► (Id Se , - Cj-,)- 1 G AA(S ex ) : 

this completes the proof. ■ 



7 Model RH Problem 

In this section, we reduce the evaluation of the integrals in Lemma 6.3 to three RH problems 
on R which can be solved explicitly. 
For k G {A,B,C}, define 

mk p (x) t j f m^-c^im^iO dt 

m iX) - 1 + k (F^ 2^' Xe£ ^ k - 

Now, from Theorem 3.1, we find that m k °(\) solves the following RH problem, 

m k +{X) = m k °{X)(I-wf)- 1 (I + wf), A G E fc , 
m fe °(A) = I- + <D{\- 2 ) as A ^ oo, AGC\£ fc . 

Substituting into (79) of Lemma 6.3 inequalities (85) and (86) (and their analogues for T,a 
and Sc), we obtain 

^. t )=7^{(^)- 2 (<°)2i+(^)- 2 (^r)2i+^)- 2 (mr)2i}+o(^). (9i) 

We consider in detail only case B. Let us introduce the function 

V{\) = m B °(A)(0(A))- 1 e- i # CT3 , A G C \ S e:c , 
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where (j>(X) is defined in (87), and notice that it is holomorphic in C \ R; in particular, it 
has no jumps across T, ex \ R. Across R, oriented from — oo to +00, using (87)-(90), one 
finds that T>(X) solves the following RH problem: 

V+(X) = P_(AMAo), A G R, (92) 
V(X) = {(l-^f- ) + 0(X- 2 )){2\y^e-^ as A -> 00, A G C \ R. (93) 

Since, according to (92) and (93), det(X>(A)) is a holomorphic, bounded in C function, 
it is, by Liouville's theorem, a constant. From (92), one finds that, det(P(A)) = 1. By 
differentiating (93), we see that d\V(X) also solves (92), so that d\V(X)(V(X))^ 1 is an 
entire function of A. Using (93), 

9 A P(A)(P(A))- 1 = -iXa 3 - i[a 3 , mf°] + C(A _1 ). 

Applying Liouville's theorem to the left-hand side of the last equation, we arrive at the 
following ODE for V(X), 

(d x V + iXa 3 V) = (3V, (94) 

where 

(3 = -i[a 3 ,mf °] = 02\o- + /3i20+; 

hence, 

( m f°) 21 = -ip 21 . 

It is well known || 26, 27 1 that the fundamental solution of (94) can be written in terms 
of the parabolic-cylinder function, D a {-) ]p4| ; in particular, 

£> + n(A) = 2-e" — A^v^Ae" — ), P +22 (A) = 2~t e~ D_ iv {y/2Xe~^), 

V +12 (X) = ((3 2 i)- 1 2-^e 1 r{d x D_ iv {V2Xe- i f) - iXD^i^Xe^ )}, (95) 
©+21 (A) = ((3 12 r 1 2^e-^{d x D lu (V2Xe~^)+iXD iu (V2Xe-^)}, 

and 



P_u(A) = 2-e-D iu (V2Xe-), D_ 22 (A) = rTe-TD^^AeT), 
Z?-12(A) = (Ai)-^"^ e-^{a AJ D^(V2Ae¥) - iXD^V^Xe^)}, (96) 
©-21(A) = (/3i 2 )- 1 2^e^{a A A,(^2Aex) + iAA,(V2AeX )}, 

where /?i 2 /3 2 i = 2za Now, substituting (95) and (96) into (92), we find that 



-Pur(Xo) = 2 i »e--W(D iu (V2Xe-),D lu (V2Xe- — )), 

where W(f,g) (= /(A)^ - g(X)^-) is the Wronskian of / = D iv {^/2\e*) an d 5 
A*(V2Ae-¥), 

W(A,(V2Aef ), A,(V2Ae-¥)) = V§|, 



and r(-) is the gamma function [24]; hence, 



B o 2 w ^e 2 e 4 
(mi )21 = r^)r(Ao) • (97) 
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From the o~ 3 symmetry reduction for m(A), i.e., m{— A) = o"3m(A)o"3, we have that, 

(mf ) 2 i = (mfV (98) 

In the case C, r(0) = r(iO) = 0, so that the corresponding model RH problem (92)-(93) 
(Ao = 0) has identity conjugation matrix, v(0) = I, whence, 

(mfhi = 0. (99) 

Hence, substituting (71), (97), (98), and (99) into (91), and recalling that Q(x,t) = P(x,t), 
we obtain (16), (18), (19), and (21) in Theorem 3.2. 



8 Asymptotic Evaluation of m(x, t] 0) 

In order to complete the proofs for q(x, t) and u(x, t) given in Theorems 3.2 and 3.3, we 
need to evaluate m(x, t;0) asymptotically (Lemma 2.1, (11)). 

Proposition 8.1 IfQ(x,Q) G S(R), then m(x,t;0) = exp{-%i \Q(£,t)\ 2 d£} . 

Proof. From the proof of Proposition 2.2, it follows that, m(x,t;0) = *$>(x,t;0) = 
exp{-^J^ o \Q(£,t)\ 2 d£}, for some x G R. Since Q(x,0) G S(R), it follows from Defini- 
tion 2.1, Proposition 2.5, and the definitions of m(x,t; A) for ^(A 2 )^ given in Lemma 2.1 

that, lim m(x,t;0) = I: note also that, m + (x,t;0) = m_(x,t;0) = m(x,t;0), since 

x— >+oo 

G(0) = I (Lemma 2.1); hence, exp{-^ f+™ \Q(£,t)\ 2 d£} =1. ■ 
Proposition 8.2 

(WQWc^R)) 2 = i (io +QO ln(1+l ; M|2) ^ - io +o ° ln(1H ; M|2) ^) . 



Proof. For Q(x,0) G 5(M), from Definition 2.1, Proposition 2.5, and the definitions of 
m(x,t;A) for 9(A 2 )>0 given in Lemma 2.1, one shows that, lim m(x,t;0) = (a (O))' 73 , 

^ X — > — oo 

where a~(0) = expjjp ln ( 1 + r M r — ^he p rQ of now follows from Proposition 2.5 
(r + (fi) = —r~(jj)), Lemma 2.1 (r(p) = r~(/j,)), and Proposition 8.1. ■ 

Lemma 8.1 As t^+oo and x^—oo such that Ao = \\J ~ f > M, 

((m-Hx,t;0)) n ) 2 = exp{f Uo° d» - / +o ° d»)} + o(*$). 

Proof. Writing |Q(£,t)| 2 ^ = -(| \Q\ | £2(R) ) 2 + J*^ |Q(£,t)| 2 ^, the proof follows 
from the asymptotic expansion of Q(x,t) given in Theorem 3.2, (16), (18), (19), and (21), 
the fact that r(A) G S(T), and Proposition 8.2. ■ 

Remark 8.1. It is possible to prove the estimate in Lemma 8.1, £?f^|^, without 
reference to the conserved quantity in Proposition 8.2 by using only the asymptotic results 
for Q(x,t): this fact was first pointed out by Ablowitz and Segur Q in connection with 
their studies of the NLSE. To follow this paradigm, one needs the asymptotic expansion of 
Q(x,t) when Ao — 0; however, we did not prove the aforementioned asymptotics for Q(x,t) 
(Remark 3.3). 

Corollary 8.1 As t — > +oo and x — > —oo such that Ao = \ \J ~ f > M , the solution of the 
DNLSE is given in Theorem 3.2, Eqs. (17)-(19) and (22). 
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Proof. Consequence of (11) in Lemma 2.1 and Lemma 8.1. ■ 

Corollary 8.2 As t^+oo and x^+oo such that A = \J\{j — \) >M an d j>\> s<ElR>o, 
the solution of the MNLSE is given in Theorem 3.3, Eqs. (32) and (33). 

Proof. Follows from Proposition 2.4 and Corollary 8.1. ■ 
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